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Abstract. We study photon-ALP conversion by resonance effects in the magnetized plasma 
of galaxy clusters and compare the predicted distortion of the cosmic microwave background 
spectrum in the direction of such objects to measurements of the thermal Sunyaev-Zeldovich 
effect. Using galaxy cluster models based on current knowledge, we obtain upper limits on 
the photon-ALP coupling constant 5 of < 0(10“^^ GeV“^). The constraints apply to the 
mass range of 2 • 10“^^ eV < utalp ^ 3 • 10“^^ eV in which resonant photon-ALP conversions 
can occur. These limits are slightly stronger than current limits, and furthermore provide an 
independent constraint. We find that a next generation PRISM-like experiment would allow 
limits down to S' ~ 0(10“^^ GeV“^), two orders of magnitude stronger than the currently 
strongest limits in this mass range. 


Keywords: axions, cosmic microwave background 


Contents 


1 Introduction 1 

2 Framework of photon-ALP oscillations 2 

2.1 Resonant photon-ALP conversion and its effect on the CMB temperature 2 

2.2 Comparing the conversion probability with the tSZ-parameter 5 

2.3 Galaxy cluster model and multiple resonances 7 

3 Results 9 

3.1 Coupling constant constraints 9 

3.2 Future perspective 12 

4 Discussion and Conclusions 13 

A Conversion probability in case of multiple level crossings 14 


1 Introduction 

While a lot of effort is put into the search for new physics at large accelerators like the LHC 
at CERN, another approach is to instead search for new physics at very low energy scales 
and small couplings. In this context, the axion is one of the best known candidates: it was 
introduced in 1977 by Peccei and Quinn to solve the strong CP problem [1], but despite 
all effort, it has not been found yet, and both its mass and its coupling to photons are 
still unknown. In addition to the axion, several extensions of the standard model predict 
similar particles, so called “axion-like particles” (ALPs) (see, e.g., [2] for a review). But 
while the axion must satisfy a certain relation between its mass and its coupling to photons 
in order to solve the strong CP problem, in general there is no relation between mass and 
coupling constants for ALPs. These ALPs have been suggested to explain several physical 
phenomena [3], such as the anomalous gamma-ray transparency [4, 5], the soft X-ray excess 
from the Coma cluster [6], or dark matter [7-10]. Interestingly, both the anomalous gamma- 
ray transparency and the soft X-ray excess from the Coma cluster may be explained with 
ALPs in a similar parameter region: the gamma-ray transparency can be resolved with a 
photon-ALP coupling constant g > GeV“^ and an ALP mass rritp < 10“^ eV [11], 

while the soft X-ray excess can be explained with g > 10“^^ GeV“^ and < 10“^^ eV [12]. 
In this work, we propose a method to improve current limits in the mass region of 10“^^ eV 
< < 10“^^ eV, therefore reducing the parameter space available for explaining the soft 

X-ray excess and the gamma-ray transparency with ALPs: for a suitable ALP mass, CMB 
photons crossing a galaxy cluster can undergo resonant photon-ALP conversion inside the 
cluster’s magnetic field, therefore distorting the black-body spectrum of the CMB. Galaxy 
clusters are one of the few large scale astrophysical objects with known magnetic helds which 
allows to derive constraints not only on the combination gB from distortions of the CMB [13] , 
but also on g itself. In particular, we use observations of the thermal Sunyaev-Zeldovich effect 
by OVRA, WMAP, MITO, and the Planck satellite and use these to obtain limits on the 
coupling between ALPs and photons. While current data leads to limits only slightly better 
than the ones obtained from SN1987A, a future PRISM-like experiment might significantly 
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improve current limits. The currently strongest limits in this mass-region are derived from 
the absence of a gamma-ray flash at the time of the SN1987A and limits the coupling to 
photons to g < 5.3 • 10“^^ GeV“^ [14]. 

This paper is structured as follows: in section 2 we introduce the framework of resonant 
photon-ALP conversion inside galaxy clusters and describe the cluster models used. The 
resulting constraints are presented in section 3. This section also contains an estimate of the 
expected sensitivity of a PRISM-like experiment. In section 4 we discuss the results and how 
they may be further strengthened. In the appendix we present in detail how the multiple 
level crossing has been calculated. 

Throughout the paper, we set c = h = ks = 1- We denote spatial vectors with bold face 
symbols. 


2 Framework of photon-ALP oscillations 

In this chapter, we derive an expression for the conversion probability from photons to 
ALPs and compare the corresponding temperature change with observations of the ther¬ 
mal Sunyaev-Zeldovich effect. Using suitable models for the profile of galaxy clusters, we 
then determine upper limits for the coupling constant between photons and ALPs. 


2.1 Resonant photon-ALP conversion and its effect on the CMB temperature 

For the deduction of the conversion probability we closely follow [13]. Axion-like particles 
(ALPs), in this work denoted by cj), are pseudoscalar bosonic particles, that couple to photons 
through the interaction Lagrangian [15] 

£ = = ( 2 . 1 ) 


where Fa/ 3/2 is the dual of the electromagnetic field strength tensor, B and E 

are the magnetic and electric field, respectively, and g denotes the axion-photon coupling 
constant. 

In an external magnetic field, this interaction Lagrangian is well known to produce effective 
mass-mixing between photons and ALPs. The new propagation eigenstates are then rotated 
with respect to the interaction eigenstates by an angle 6 given by [16] 


sin 20 


2gBoj 


cos 20 = ^ =. 

+ {2gBujY 


( 2 . 2 ) 


Here, uj denotes the energy of the photon, B is the component of the magnetic field per¬ 
pendicular to the propagation direction of the photons, is the ALP mass and g is the 
coupling constant as before. From here on, we will refer to 0 as magnetic mixing angle. 
Using typical values for the parameters in our study, the relevant dimensionless parameter 
in these expressions reads 


~ 1.38 • 10-6 3 B uj / 10 eV \ 

10-12 GeV ^ gG 10 ^ eV V J 


(2.3) 


For the parameter ranges considered here this will always be much smaller than unity. 
The misalignment between interaction eigenstates and propagation eigenstates will produce 
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photon-ALP oscillations with a wavenumber given by [16] 


k = 


2io 


(2.4) 


Inside a plasma, photon-ALP mixing will be modified: the refractive properties of the plasma 
lead to a non-trivial dispersion relation, which can be parametrized by an effective photon 
mass m^. In this case, for a given magnetic field the effective mixing angle in the plasma 9 
is related to the mixing angle 9 at vanishing charge density by [17] 


where ^ is defined as 


sin 20 


cos 29 


sin 20 

(2.5) 

[sin^ 20 -|- (cos 20 — ^)^]^/^ ’ 

cos 20 — ^ 

( 2 . 6 ) 

[sin^ 20 -|- (cos 20 — ^) 2 ]i /2 ’ 

/ \ 2 

= cos 20 ^ . 

(2.7) 



If in some region in space the resonance condition 


= 7710 


( 2 . 8 ) 


is satisfied, one has 0 —?■ vr/d and resonant photon-ALP conversion is possible. 

In this study, we will consider such resonant photon-ALP conversion occnrring in galaxy 
clusters, where one has both external magnetic fields and a nonzero free electron density. 
Due to the high ionization fraction in the intra-cluster medium, contributions from the scat¬ 
tering off neutral atoms can be neglected and the effective photon mass is given by the plasma 
frequency [18] 


2 2 ^TTO 

m ~ CUp = -Ue, 

' rUe 


(2.9) 


where a is the fine structure constant, nie is the electron mass and Ue is the free electron 
density. Using this relation, one can rewrite the resonance condition (2.8) to 


m 0 = 3.72 • 10"^^ eV. 


( 2 . 10 ) 


As the free electron density cannot reach arbitrary values inside a galaxy cluster, the reso¬ 
nance condition can only be satisfied for a certain range of ALP masses m^. Assuming for 
example a minimal density of 0.3 m“^ and a maximal density of lO'^ m“^, resonant photon- 
ALP conversion will only occur for 2 • 10“^^ eV < < 3.7 • 10“^^ eV. The first value is 

the average baryon density in the universe at redshift zero and the second value is a typical 
density in the core of a galaxy cluster. 

The distance between the photon production and the resonance as well as the distance be¬ 
tween the resonance and the detection is much larger than the oscillation length causing an 
incoherent superposition of the oscillation patterns. In this case, the transition probability 
is given by [19] 





COS 200 cos 29 D, 


( 2 . 11 ) 
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where 6q is the effective mixing angle at production, 9d is the effective mixing angle at detec¬ 
tion and p is the level crossing probability. Therefore, a transition from a medium dominated 
to the vacuum state corresponds to cos20o — —1 and cos 29 d — cos 20 or vice versa and for 
0 <C 1 one obtains a conversion probability Py^fj, close to unity for p <C 1, corresponding to 
an adiabatic transition. 

As it was argued in [13], the high plasma density at the time of the creation of the CMB 
photons leads to a value of 9i close to 7r/2. Using typical values of the free electron density 
in the solar system, one can see that also 0z) is very close to 7r/2. A more detailed discussion 
of this point can be found in the appendix. 

For our range of ALP masses, there is not only one, but several resonances: the first reso¬ 
nance occurs when the free electron density decreases due to the cosmic expansion. Inside 
the transversed cluster, there are several resonances: additionally to the two resonances due 
to the increasing (decreasing) electron density when entering (leaving) the cluster, density 
fluctuations enable even more resonances. Finally, there are also resonances when the pho¬ 
tons enter the Milky Way. 

The level crossing probability pi for a single resonance i is given by the Landau-Zener ex¬ 
pression [20] 

Pi ~ exp(—27ri?A:sin^ 0res), ( 2 - 12 ) 

where k is the oscillation wavenumber given in eq. (2.4), 0res is the magnetic mixing angle 
given by eq. (2.2) at this resonance and R is the scale parameter defined as 


R = 


din 

dt 


-1 


t=t 


res 


(2.13) 


The level crossing probability pi takes into account the deviation from adiabaticity of the 
photon-ALP conversion in the resonance region. One has p* ~ 0 for a completely adiabatic 
transition and Pi = 1 for an extremely non-adiabatic one. 

The Landau-Zener expression (2.12) only holds for the case when the free electron density 
varies linearly during the resonance. For 0res the resonance half-width is, according to equa¬ 
tion (2.5), 6^ ~ sin20i.es, corresponding to a resonance width in the density scale of 


An^ ~ sin(20res) — n 


j^2gBuj 

9 

mi 


~ 10"®n^, 


(2.14) 


where is the resonance density defined by eq. (2.8). Due to this extremely narrow 
resonance region, the approximation of a linear density change is very well fulfilled. 
Expanding the sine in eq. (2.12) and approximating k ~ m'^j{2uj), the exponent becomes 

27ri2A:sin2(0,es) ^ 0(10-®), (2.15) 

ml 


for typical values used in this work. As derived in Appendix A, for such small exponents, 
the total level crossing probability becomes 


p=n 

i 


Pi ^ 


^ g^BjRjUjTT 

^ ml 

i 0 
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where the sum is over all resonances i. Together with the obtained expressions for cos20o, 
and cos20 £), the conversion probability then is 




g^BfRiOjTT 




m. 


g^UlTT 

2 

ml 




(2.17) 


In the present study we will neglect the first resonance because of the unknown magnetic 
field, and the resonances in the Milky Way because of the small scale parameter R. These 
approximations are conservative, since they tend to underestimate the actual conversion 
probability, as it will be discussed in Appendix A. 

The conversion of photons into ALPs will always reduce the number of photons in the beam, 
cansing an apparent temperatnre decrease. The intensity for a given photon energy and 
temperature I{uj,T) can be related with the apparent temperature change over 


_ dI{uj,T) ^ AI 
~ dr “ ^ 


—^ AT{u]) ~ 


AI{u,T) 

P{u^,T) ■ 


(2.18) 


The intensity of the photon beam after crossing the galaxy cluster is reduced by a factor 
of 1 - giving AI{uj,T) = -P^^^Io{uj,T). With I{uj,T) oc [exp(t<;/rcMB) - 1]”^ one 


arrives at 


AP<y —^^(cu) P'y — 


T'i 


CMB 

UJ 


1 — exp — 


Tcmb ) 


(2.19) 


for the apparent temperature change in dependence of the photon energy due to resonant 
photon-ALP conversion. 


2.2 Comparing the conversion probability with the tSZ-parameter 

From the Planck 2015 data [21], we have information about the temperature differences in the 
directions of galaxy clnsters. These temperature differences depend on the freqnency observed 
and are related to the (by definition freqnency-independent) thermal Sunyaev-Zeldovich (tSZ) 
Compton parameter y by the expression [22] 

Artsz(w) = /(w) Tcmb y, ( 2 . 20 ) 


with /(w) = (w/Tcmb) • coth[c<;/(2rcMB)] — 4. The function /(w) is negative for values 
of u smaller than 3.83 TcmB) otherwise the function is positive. In figure (1) the relative 
temperature change due to the thermal SZ effect as well as due to photon-ALP conversion 
is shown. As mentioned above, the relative temperature change dne to the photon-ALP 
conversion is always negative, because the effect always removes photons from the beam. 
The tSZ-effect, in contrast, creates a negative temperature change for small energies (w < 
3.83 Tcmb) and a positive temperature change at higher photon energies. 

For the Coma cluster, detailed measurements of its thermal SZ-effect for photon energies in 
the range of 10““^ eV < w < 1.1 • 10“^ eV exist [23]. These values are presented in table (1). 
In this case, one can build the reduced y^-function 


x^{y,g) 


1 


N r 


N -2 


E 




(J, 


exp 


( 2 . 21 ) 


where are the observed temperature changes at photon energy Ui, Oi are their standard 

errors and 

^j.theo(y^^) ^ AT^^^{uJi,g) -h ATtsz(wi,y), (2.22) 















Figure 1. Comparison of the spectral profile of the tSZ effect (with y = 10“®) and the photon-axion 
conversion. The a;—axis shows the photon energy, the y—axis shows the relative distortion of the 
effective CMB temperature, AT/Tcmb in units of 10“®. For this plot, the parameters g = 5 ■ 10“^® 
GeV“^, B = 2 fiG, R = 0.5 Mpc, = 10“^® eV were used. Different values for these parameters 
change the normalization of the effect but not its dependence on the photon frequency. The filled 
circles refer to the centers of the frequency bands of the Planck-mission. 


Experiment 

w [10“^ eV] 

A Ttsz [ pK ] 

OVRA 

1.32 

-520 ± 83 

WMAP 

2.51 

-240 ± 180 

WMAP 

3.87 

-340 ± 180 

MITO 

5.91 

-184 ± 39 

MITO 

8.85 

-32 ± 79 

MITO 

11.25 

172 ± 36 


Table 1. Measurements of the thermal SZ-effect in the Coma cluster. 

is the prediction by the theory, with and ATtsz given by eq. (2.19) and eq. (2.20), 

respectively. For each ALP mass rricj), we calculate the x^-function in the (y, y)-parameter 
space and determine the limit on g as the largest value of it still inside the respective confi¬ 
dence interval. 

Often, e.g. in the case of the Planck-mission [21], only the y-parameters of galaxy clusters 
are easily available, but not the temperature changes due to the thermal SZ-effect at different 
frequencies. In such cases, a simple and conservative bound can be obtained when assuming 
that the magnitude of the temperature change due to the photon-ALP conversion must be 
smaller than the magnitude of the temperature change due to the tSZ-effect: 

< lATtszl (2.23) 

Using equations (2.17), (2.19), (2.20), and solving for the coupling constant g, one arrives at 

/ \ 1/2 
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where we defined 


h{u}) = 


4 — X • coth(x/2) 


1/2 


LO 


X = 


Tcmb 


(2.25) 


rcMB7r[l - exp(-x)] 

Note that this approach does not take any non-resonant conversion into account, and therefore 
conservatively overestimates g. 


2.3 Galaxy cluster model and multiple resonances 

Both the magnetic field strength as well the scale parameter at resonance depend on the 
free electron density inside the considered galaxy cluster. On top of the smooth, large- 
scale electron density profile, galaxy clusters exhibit smaller turbulent contributions in the 
electron density. For example, in the Coma cluster, density fluctuations of ~5% on scales of 
~30 kpc and ~(7-10)% on scales of ~500 kpc have been observed [24]. As estimated before, in 
terms of the electron density variation the resonance is extremely narrow, Aue/ures ~ 10“®, 
such that small turbulent contributions can already cause several distinct resonances. We 
therefore have to evaluate the expression from eq. (2.17), where the index i labels 

the individual resonances. 

We assume that the magnetic field strength follows the free electron density, such that 


(|B(r)|) ~ Bmax[ne(r)/nmax]'', (2.26) 

where g < 0(1) and Bmax will be specified later. Due to the narrowness of the resonance 
region, one can approximate the magnetic field as being constant during the resonance. Fur¬ 
thermore, the magnetic field during the resonance is completely determined by the resonance 
density nj-^s, such that all resonances will experience the same field strength. Averaging over 
several resonances, one therefore obtains 

( 2 . 27 ) 

i i 

where the factor of 2/3 accounts for the fact that only the two transverse components of the 
magnetic field enter the conversion probability (2.17). 

To investigate the consequences of turbulent density contributions, we will consider an elec¬ 
tron density with a dominant smooth component rig and a spectrum of modulations with 
wavenumber k, amplitude 5{k), and phase c/fc, such that 

ne(r) = ns(r) • ]^[1 -b 5(/c) sm{kr + /)*.)]. (2.28) 

k 

We will assume some upper cutoff /cmax induced by viscous damping and 6{k) 1 as indicated 

by observation. As Auj-es/nres ^ 1, one can neglect “incomplete” resonances (where the 
density has an extremum), but always assume the electron density to vary linearly during 
the resonance. The scale parameter Ri, defined by eq. (2.13), of an individual resonance i 
at radius r* then is 

^ dln[ng(r)] 5ik)k cos{kr + ^k) , 

dr ^ ^ l+S{k)sm{kr + cl)k) r=n ^ ’ 

Using <5(500 kpc) ~ (7 — 10)%, <5(30 kpc) ~ 5% from the Coma cluster, and assuming a 
power law 6{k) oc k~^ provides ^ ~ 0.12...0.25, implying 6{k)k oc /jO'75...o.88^ scale 
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parameter will therefore be dominated by the contribution from the largest wavenumber not 
affected by damping, while, due to the random phases the contributions from larger 
scales approximately average out and can be neglected. With 6{k) <C 1, and denoting the 
dominating wavenumber as kdom, one thus has 


Ri ~ 


(iln[ns(r)] 


dr 


+ (5(A:(joni)/C(jom COs(fe(^om^i “1“ 4^k) 


-1 


(2.30) 


The spatial width Ar containing all the resonances can be estimated by the width, within 
which the smooth profile varies by a factor of (1 TJmax), where (5 max is the largest modulation 
amplitude. Explicitly, Ar ~ 2<5max^res/|^sl) where = dns{r)/dr\j-es is the derivative of the 
smooth profile, evaluated at resonance. As ^(/cdom) < <5(A:max), one can conservatively set 


Ar ~ 2(5(/Cdom)^res/|^sl) 


(2.31) 


neglecting a factor of order unity. Including additional modes would enable resonances within 
an even larger region, making this statement only more conservative. 

The number of resonances can then be estimated by 

n^domAr 2<5(A;dom)^dom^res 

27r 7r|^sl 


(2.32) 


where the factor of two in the hrst equality arises because there are two resonances per 
complete period. Note that there must always be at least one resonance (if n^es < W ma x) for 
radially incoming (outgoing) photons due to the increasing (decreasing) electron density. 

As N has to be a natural number, this formula is only a good approximation for N ^ 1, while, 
for A « 1, we expect an uncertainty of order unity. The assumption N ^ 1 implies certain 
conditions on the density fluctuations: density fluctuations of ~ 5% on scales of 30 kpc have 
been detected in the Coma cluster, providing N > 2, while projection effects preclude strong 
limits for smaller scales [24]. Rotation measures, however, indicate that the Coma cluster 
contains magnetic fields with coherence lengths down to ~ 2 kpc [25]. Due to turbulence, 
one expects density fluctuations on similar length scales. Again using 6{k) oc k~^ with 
^ ~ 0.12...0.25, one obtains 6{2 kpc) ~ (2.5...3.6)%, while N I implies 6{2 kpc) 0.2%. 
In the last inequality, we assumed the /3-model introduced below as a smooth profile and 
the parameters of the Coma cluster. In the same cluster, and <5(2 kpc) ~ 3%, one obtains 
N > 15. For such a high number of resonances, the second term in (2.30) dominates, and 
one can approximate the sum over all resonances by averaging the trigonometric functions, 
leading to 

^ ~ ^ |(5(/Cdom)fcdom COs(/CdomD + </>fc) 

i i 


~ N 


vr 




= (2.33) 


where (j cos(x)j) = 2/tt and eq. (2.32) for N has been used. Rs is the scale parameter one 
obtains from the smooth profile without density modulation. 

If N is 0(1), the scale parameter (2.29) becomes Ri ~ Rs/N. One could then still approxi¬ 
mate Ri ~ Rs, inducing a relative error of < 0(1). 

As we have seen, the exact dominating wavenumber fedom and corresponding amplitude 
'5(fcdom) do not influence the transition probability as long as N ^ 1. We will therefore 
not specify them in any more detail and approximate 

^^BiRi ~ (2.34) 
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Independent of the multiple resonances due to the turbulent structure, there is one region of 
resonance when the photons enter the galaxy cluster, and one region of resonances when the 
photons leave the cluster. Thus, an additional factor of 2 has to be included in the conversion 
probability. In total, one therefore obtains 


-P7—— 




(2.35) 


For numerical calculations, we will consider a /3-profile as the dominant smooth profile 

ns(r) = max nmax[l + , (2.36) 

where (3 is 0(1), Umax is the free electron density in the cluster center, and rc is the core 
radius of the cluster. We include the average cosmological electron density no as a lower 
boundary for the electron density, and therefore, due to eq. ( 2 . 8 ), a lower cutoff for the 
ALP-masses able to undergo resonance. 

We furthermore will focus on two clusters: the Coma cluster and the Hydra A cluster. In 
[25], rotation measure images have been used to determine the Coma clusters magnetic field 
strength as well as the parameter rj, defined by eq. (2.26). In this analysis, degeneracy 
between Hmax and rj has been found; a larger Hmax implies a larger g and vice versa. The 
best fit gave Hmax = 4.7/rG, t] = 0.5, while Hmax = S.dgG, rj = 0.4 and Hmax = 5.4/iG, 
r] = 0.7 are still within 1 a. To illustrate the dependence of our approach on Hmax and rj, we 
will perform our analysis with these three pairs of values. We furthermore adopt the values 
Tc = (291 ± 17) kpc, Umax = (3.44 ±0.04) • 10^ m“^, and /3 = 0.75 ±0.03 from the same work. 
In contrast, the Hydra A cluster is a cool-core cluster and exhibits magnetic field strengths 
of 6 jiG coherent on scales of 100 kpc and magnetic field strengths of 30 jiG on scales of 4 
kpc [26] . The electron density in the cluster center is Umax — 10^ m“^ and the core radius is 
Tc ~ 130 kpc [27]. We also adopt the value (3 = 1 from [27], while different values of rj have 
been used in the same work: mostly, g = 0.9 has been used, but also g = 1/2, g = 2/3, and 
g = 1 have been considered. Due to its strong influence on the possible limits on g, we will 
work with g = 0.9 as well as with g = 2/3. The former value is suggested by observations 
of Abell 119 [28], while the latter is predicted by flux conservation and is closer to the value 
observed in the Coma cluster. 

For both clusters, we will assume y ~ 10“®, a typical value of the Compton parameter 
observed in galaxy clusters [ 21 ]. 


3 Results 

3.1 Coupling constant constraints 

In figure (2), the limits obtained from the y^-analysis of the Coma cluster are presented. Due 
to the strong influence of g on the magnetic field strength in the outer regions of the cluster, 
the limits from three different pairs of values for (Hmax; v) shown. These three pairs of 
values are the ones already presented in the previous section: (4.7/iG, 0.5) provides the best 
fit to the observed rotation measures, while (3.9//G, 0.4) and (5.4/rG, 0.7) are the lower and 
upper limits, respectively, at la confidence level. 

For the best-fit model from [25], the obtained limits are up to a factor of 5 stronger than 
the limits derived from SN1987A. For Hmax = 3.9/rG, g = 0.4, the obtained limits are even 
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Figure 2. Obtained limits for different pairs of values for i?max)?7 obtained by the x^"test with 
the data from the Coma cluster. The black solid line shows the limits at 95% C.L. when assuming 
the best-fit model from [25]. Using the upper and lower limits on Umax)?? at Icr, one obtains the red 
dashed-dotted line and the green dashed line, respectively. The limits at 99% C.L. for the different 
values are indicated by the dotted lines. While the limits on the coupling constant g for high ALP 
masses (corresponding to resonances near the center of the cluster) are quite similar, the difference 
increases for smaller ALP masses (corresponding to resonances in the outer regions of the cluster). 
The shaded area is excluded by limits obtained from SN1987A. 


stronger for small ALP masses, while i?max = 4.7/iG, rj = 0.7 produces limits slightly weaker 
than the ones from SN1987A. 

In this context, a warning is necessary: the /3-model of the free electron density as used in 
[25] and adapted here is based on measurements of the Coma cluster’s X-ray emission [29]. In 
this work, the largest distance from the cluster center probed is ~15 Vc, as the signal becomes 
undetectable in the noise for larger distances. This maximal tested radius corresponds to 
Ue — 8 m“^ and ~ 1.0 • 10“^^ eV. Although the density obviously has to decrease to 
the average cosmological density, the exact prohle is not known. We will assume that the 
/3-profile holds down to the average cosmological density, but one should keep in mind that 
for < 10“^^ eV, the limits are obtained under this assumption. 

In figure (3), a comparison of the limits obtained by the x^-test (2.21) and the limits ob¬ 
tained by the temperature comparison (2.23) is presented. For the limits by the temperature 
comparison, the value of g was determined according to eq. (2.34). Additionally, we included 
Gaussian error propagation to estimate the uncertainty and to obtain limits at different 
confidence levels. We adopted the uncertainties given by [25], i.e. = 16%, 

'^n.max/^max = 1-1%, CTp/P = 4%, Ur^lvc = 6 %, as well as (Toj/o; = 20 %(low frequency instru¬ 
ment)/33% (high frequency instrument) from [30]. Although the total conversion probability 
(2.34) does not dependent explicitly on N anymore, the multiple resonances induce an ad¬ 
ditional uncertainty. When assuming X > 40 (see above), and oc y/N, one obtains 
(Jn/N ~ 16%. As Yli Ri oc N (X Rp oc Tc, we absorb this uncertainty into Or^jrc and conser¬ 
vatively set Or^jrc — 25%. Finally, we conservatively set Oyjy = 50%. The error is usually 
dominated by the uncertainty of the photon energy cj; only for very low and for very high 
photon energies, the uncertainty is dominated by Oy. 
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Figure 3. Comparison of the limits at 95% C.L. obtained by the y^-test (2.21) and the direct 
temperature comparison (2.23) with different photon energies. For the photon energy oj = 8.97-10“^ 
eV (green solid line), the function h{ijj) in eq. (2.25) becomes almost zero. But the rapid change of 
the function causes a high uncertainty, effectively even weakening the obtained limits. For all other 
(center-) frequencies measured by Planck, the obtained limits lie inside the red/light gray area, where 
w = 35.44 • 10“"^ eV gives the weakest bounds and uj = 5.91 • 10“"^ eV gives the strongest bounds. 
These bounds are approximately 20% weaker than the limits obtained by the x^-test (black solid 
line), proving that the direct temperature comparison provides realistic results. The blue/dark gray 
shaded area is excluded by limits obtained from SN1987A. For this plot, the data of the Coma cluster 
and ?7 = 1/2 has been used. 


In order to avoid the singular behavior of for n^es ric, we exclude the innermost region 
with Ures/nc > 10%, and, due to its large influence, we fix r? = 1/2. This plot should therefore 
demonstrate the robustness of the obtained limits with respect to astrophysical uncertainties. 
The used photon energy determines the numerical value of the function h{uj), see eq. (2.25). 
This function is usually of order unity, but reaches zero for ui = 3.83 Tcmb- Naively, one could 
expect to obtain arbitrarily strong limits when simply using a photon energy close to this 
value. This is, however, unphysical: the Planck high-frequency channels have bandwidths of 
Azz/zz ~ 0.33 [30], meaning that every frequency map is actually an average over a range of 
frequencies. One therefore also would have to take an appropriate average over the function 
h{u)), preventing arbitrarily small limits. 

A comparison of the limits obtained from the Coma cluster and from the Hydra A cluster 
using eq. (2.34) is shown in figure (4). For the Hydra A cluster, we have assumed 25% uncer¬ 
tainties for all cluster parameters, i.e. rc, Umax) H ma x, /3, and kept ay = 50% as before. The 
higher central electron density in the Hydra A cluster allows higher ALP masses to undergo 
resonant conversion, therefore slightly expanding the mass range accessible for the method 
presented here. The higher value of r] leads to a faster decrease of the magnetic field strength 
with increasing radius, such that the obtained limits become weaker for smaller ALP masses. 
For illustration, we also show the limits for the Hydra A cluster with r/ = 2/3 and the Coma 
cluster with rj = 0.7: in this case, the radial decrease of the magnetic field is very similar, 
while the higher value of the central magnetic field in the Hydra A cluster leads to a higher 
conversion probability and therefore slightly stronger limits. 
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Figure 4. Comparison of the obtained limits at 95% C.L. for two different clusters: the Coma 
cluster and the Hydra A cluster. Additionally to each labeled line, the limits at 99% C.L. are shown 
as dotted lines. The higher central density of the Hydra A cluster allows limits for slightly higher 
masses, while the strength of the limits strongly depends on the assumed values of rj. For this plot, 
w = 5.9 • 10““^ eV has been used. 


3.2 Future perspective 

In a more detailed study, one would have to simultaneously fit several contributions to the 
recorded data, e.g. the tSZ effect, thermal dust and synchrotron radiation. Including photon- 
ALP conversion in this procedure, one would then obtain limits on the coupling constant g. 
To estimate the possible limits with this approach, we restrict ourselves to a simpler ap¬ 
proach: we simulate a tSZ-signal according to eq. (2.20), where we use the uncertainties of 
the Planck-experiment, multiplied with a factor referred to as “error penalty”. This error 
penalty parametrizes the additional uncertainty induced by subtracting the foreground emis¬ 
sion. In a second step, we perform a x^-cinalysis, according to eq. (2.21), where we fit both 
a tSZ-signal and photon-ALP conversion to the simulated signal. We thus obtain an upper 
limit on the sensitivity for the coupling constant g. 

We use the temperature sensitivities described in [31] and [30], where the sensitivities range 
between AT/T ~ 2.2 • 10“® for v = 143 GHz and AT/T ~ 6 • 10“^ for v = 857 GHz. In figure 
(5), we show the possible limits with error penalties of 1, 5, and 10, where the parameters 
of the Goma cluster have been used and we averaged the limits obtained from ten different 
simulated realizations. 

One also can extent this approach to proposed future experiments for highly sensitive CMB 
observation like PRISM [32] or PIXIE [33]. PRISM is proposed to have 32 broad-range fre¬ 
quency channels as well as 300 narrow frequency channels covering the range from 30 GHz to 
6000 GHz. The simulated 4-year sensitivities range from = 3.6x 10“^^ Wm“^Hz“^sr“^ for 
frequencies between 30 GHz and 180 GHz to 5Iu = 1.6x10“^® Wm“^Hz“^sr~^ for frequen¬ 
cies greater than 3 THz. Performing the same analysis as before, we arrive at sensitivities of 
> orders of magnitude below the currently strongest limits for this mass 

range. This result is also displayed in figure (5). PIXIE is proposed to have 400 channels cov¬ 
ering the same frequency range as PRISM, reaching slightly worse sensitivities than PRISM. 
We therefore expect PIXIE to be sensitive to values of g very similar to the ones presented 
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Figure 5. Projected sensitivity of the presented approach when using the full Planck data or data 
from a future, PRISM-like experiment. The sensitivity of Planck with no error penalty is shown with 
the solid red line, while the dashed black line corresponds to a PRISM-like experiment. The dashed- 
dotted lines are obtained when an error penalty of 5 is used, the dotted lines for error penalties of 10. 
The shaded area is excluded by limits obtained from SN1987A. For this plot, the parameters of the 
Coma cluster have been used. 

for PRISM. 

4 Discussion and Conclnsions 

In this study, we have shown that Planck’s recent measurements of the tSZ Compton- 
parameter y [21] can be used to constrain the coupling constant g of pseudoscalar ALPs 
to photons. To this end, we compared the temperature change due to the tSZ-effect of the 
CMB photons reaching us from galaxy clusters with the temperature change due to resonant 
photon-ALP conversion. Photon-ALP conversion is most effective at resonance where it leads 
to the strongest limits on the coupling constant. On the other hand, resonant photon-ALP 
conversion is only possible for a limited range of ALP masses, typically of the order of 10“^^ 
eV; this range depends on the effective photon mass in the galaxy clusters, and therefore on 
the free electron density. 

The strength of the obtained limits depends both on the density prohle in the galaxy cluster 
as well as on the magnetic field. In our study, we used a /3-model, extended with density 
modulations, for these profiles, and typical values for the magnetic field strength, electron 
densities and the observed y-parameter for galaxy clusters. Under these assumptions, we can 
derive limits on the photon-ALP coupling constant g, which are slightly stronger than the 
existing bounds in this mass region from SN1987A [14], and furthermore provide an inde¬ 
pendent constraint. 

The scaling of the magnetic field with the electron density has a strong influence on the 
limits for smaller ALP masses. Further investigations of the magnetic field strength in the 
outer regions of galaxy clusters or the selection of suitable clusters will be necessary to solve 
this problem and might provide stronger limits. Another approach might be to consider the 
magnetized jets emitted by AGNs. 

We estimated the parameter region a PRISM-like experiment would be sensitive to and found 
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that for the mass range given by the resonance condition, sensitivities down to g > 10 “^^ 
GeV“^ are realistic. This is especially interesting, as this is part of the parameter space 
has been invoked to explain the soft X-ray excess of the Coma clnster [ 12 ] or the anomalous 
gamma-ray transparency [ 11 ]. 
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A Conversion probability in case of multiple level crossings 


In [19] it was derived that the photon-ALP conversion probability is given by 


P.y^^ = {p~^ cos( 20 o) cos(20i5). 


(A.l) 


where Oq is the effective mixing angle at the production of the CMB, is the effective mixing 
angle at the detection, and p is the total level crossing probability. When averaging over the 
distances between the individual resonances, the total level crossing probability is the sum 
of all the individual probabilities of getting an odd number of level crossings, so simply the 
classical probability results. 

We will first consider an odd number of resonances. In the case of only one resonance, the 
level crossing probability is just given by this single level crossing probability. In case of three 
level crossings, one has 


P = Pl \P2P3 + (1 -P2)(l -Ps)] + (1 -Pl) b2(l -Ps) + (1 -P2)P3] , (A.2) 

where pi, i = 1, 2, 3, are the level crossing probabilities for the individual resonances. These 
individual probabilities can be calculated with the Landau-Zener formula [19] 

Pi ~ exp(—27rRA:sin^ 0res) =: exp(— dj), (A.3) 

where R is the scale parameter defined in eq. (2.13), k is the wavenumber of the photon- 
ALP conversion (2.4) and Ores is the magnetic mixing angle given in eq. (2.2), evaluated at 
resonance. Note that this formula only holds when the electron density varies roughly linearly 
within the resonance. In our study, the exponent 5i is very small, so one can approximate 
exp(—(5j) ~ 1 — (5j. Plugging this into formula (A.2), one obtains 

3 

p = I-'^di + '^didj - 4 :-616263. (A.4) 

i=l i^j 

This result generalizes to any odd number n of resonances: 

n 

p = 1 — ^ 6i -b 0(6“^ -b higher). (A.5) 

i=l 
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In case of an odd number of resonances, one necessarily has cos(20j) < 0 and cos(20j) > 0 or 
vice versa. Parametrizing cos(20o) = =t(—1 + eo), cos(20£)) = ±(1 — eo) with |eo|, |e£)| < 1, 
one has 

cos(20o) cos( 20£)) = -1 + eo + e^) - eoeo. (A.6) 

Using the last two formulas, the conversion probability (A.l) for an odd number of resonances 
reads 

Si + (eo + cd) ^1/2 — — eoCD ^1/2 — 5^ + 0(5'^ + higher). (A.7) 

In our study, we have four resonances. In this case, the total level crossing probability 
becomes 


P = Pl{P2 b3(l -P4) + (1 -P3)p^] + (1 -P2) \P3P4 + (1 “ Ps)(1 “ P 4 )]} + 

+ (1 -pi){p2 b3P4 + (1 -P3)(l -P 4 )] + (1 -P 2 ) b3(l -P 4 ) + (1 -P3)P4]}, 


where pi are again the individual level crossing probabilities. Using formula (A.3) and <5j <C 1 
as before, one arrives at 

4 4 

p = ^ ^ + 4 • ^ {5i62S364)/6i - 8 • (5i(52(53(54. (A.9) 

2 = 1 i:^j 2 = 1 

More general, for any even number n of resonances, the total level crossing probability will 
be 

n 

p = 6i + 0(6^ + higher). (A.10) 

2=1 

In case of an even number of resonances, one either starts below the resonance density, 
and also ends below the resonance density, or one starts and ends above the resonance 
density. One can therefore parametrize cos(20o) = ±(—1 + ^o)) cos(20i5) = ±(—1 + cd) with 
bol) \^d\ < 1, and obtains 

cos(20o) cos( 20£)) = 1 - eo - Ed + eoeo- (A.11) 


Using the last two formulas, conversion probability becomes 


Si + (eo + e^i) ( 1/2 — Si 




1/2-Y,Si]+ 0{S^ + higher), (A.12) 


which is, up to terms of 0{S^ + higher), the same formula as before. 

As it was argued in [13], due to the high electron density at the time of recombination 
cos ( 200 ) = — 1 + eo is very close to —1, with 


eo — 7 • 10 


-18 


9 


B. 


cosmological ^ 


10-13 GeV 


-1 


nG 


Tcmb) 


(A.13) 


where i?cosmoiogicai is the magnetic field strength on cosmological scales and all values are 
taken today. Using a free electron density of about 1 cm-3 = 10® m-3 inside our galaxy, one 
obtains cos(20£)) = —1 + eo with 


ei) ~ 5 • 10 


-15 




10 


-14 


eV 


+ 7-10 


-12 


9 


B, 


galactic ^ 


10-13 GeV 


-1 


pG 


V 

TcmbJ 


(A.14) 
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The first term can reach values up to 10“® for the ALP-masses considered here. Therefore, 

all terms of order and (5 • e in the conversion probability can certainly be neglected and, as 
<C 1, also all terms of order 5^. The remaining result simply is 

— E + (eq + cd)/2. (A.15) 

i 

In this study, we will always make the conservative assumption that 

(A. 16) 

therefore conservatively underestimating the conversion probability, since any possible (pos¬ 
itive) contribution from cq and e^) is neglected. As we are interested in upper limits for g, 

underestimating is conservative and only makes the limits more reliable. 
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